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Abstract. A phenomenological quasiparticle model is surveyed for 2+1 quark fla- 
vors and compared with recent lattice QCD results. Emphasis is devoted to the 
effects of plasmons, plasminos and Landau damping. It is shown that thermody- 
• namic bulk quantities, known at zero chemical potential, can uniquely be mapped 

towards nonzero chemical potential by means of a thermodynamic consistency 
, condition and a stationarity condition. 

D ' 

1 Introduction 

(N ' 

Intense experimental and theoretical investigations \T\ suggest the existence of a new, deconfined 
phase of strongly interacting matter, where quarks and gluons form a fluid or gas, the quark- 
gluon plasma (QGP). If confirmed, the QGP would have existed during the Big Bang prior to 
hadronization and might be found inside of massive neutron stars. Indeed, recent results from 
the Relativistic Heavy Ion Collider (RHIC) experiments point to the formation of a quark-gluon 
O^N medium of low viscosity [2)3j4|5j . 

However, on the theoretical side, much work remains to be done. Perturbative solutions 
of QCD [6171819110111) are limited to the region of asymptotic freedom and fail for the strong 
coupling regime (e.g. in the vicinity of the pseudocritical temperature |12j of deconfinement 
T c ; at somewhat higher temperatures, say above 2T C resummation improves the convergence 
of perturbation theory noticeably [13J14J15] ). Numerical evaluations of the full theory, on the 
other hand, are still restricted to small chemical potential as being useful for the Big Bang or 
heavy-ion collisions at present RHIC top energies or future LHC energies. However, at RHIC 
bottom energies, at SPS energies and, in particular, at FAIR energies, baryon density effects 
become significant and require different approaches. 

Quasiparticle models (QPM), describing the quark-gluon plasma as assembly of essentially 
non-interacting excitations emerging from the strong interaction, have proven to represent useful 
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phenomenological parametrizations of QCD thermodynamics above T c [16jl7U8jl9l20| . At zero 
chemical potential, lattice results are described with surprising accuracy allowing the adjust- 
ment of model parameters. Thermodynamic self-consistency, supplemented by the stationarity 
of the thermodynamic potential, can then be used to extrapolate thermodynamic properties of 
systems to nonzero net baryon density. 

Our quasiparticle model (21 22 23 24 is based on the HTL approximation [25] to the 1-loop 
self- energies. This gives rise to four quasiparticle families. While quasiquarks and transversal 
gluons within the model represent excitations with quantum numbers of actual quarks and 
gluons with modified masses, quark holes (plasminos) and longitudinal gluons (plasmons) are 
quanta of collective excitations. The residues of the poles in the spectral density of the collective 
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modes vanish exponentially for momenta k ~ T, from which the main contributions to thermo- 
dynamic phase space integrals originate. Therefore, they were neglected in the previous simple 
form of the model (dubbed cQP in 26J). Additionally, damping contributions were neglected 
as they are small at zero chemical potential, fi — 0. 

The procedure of mapping the eQP results from // = into the T-/i plane is plagued 
by some ambiguities leading to non-unique solutions close to the presumed phase transition. 
Therefore, the model is restricted to sufficiently large temperatures. First attempts to include 
collective excitations into a two-flavor quasiparticle model at nonzero chemical potential have 
been made [27] and suggest that these ambiguities might vanish. In this work we show that 
both collective excitations and damping effects are necessary to preserve the self-consistency 
of the model and ensure unique solutions when extrapolating towards large baryon densities 
at moderate temperatures. In the present work, the 2+1 flavor case is considered, allowing the 
use and extrapolation of recent lattice data [28] . 

Our paper is organized as follows. Section [2] will comprise the derivation of the full HTL- 
based QPM, including collective modes and damping, as a series of approximations from QCD. 
The necessary further approximations leading to the eQP and its problems are discussed in 
section [3] The results for both models are then contrasted in section [4] and investigated in some 
detail. Finally, a conclusion is given in section [5] 

2 Derivation of the full HTL model 

2.1 The effective action 

A connection of the fundamental theory of QCD and the thermodynamic potential of the 
QGP is provided by the Luttinger-Ward formalism [29)30] as shown in [23] . Alternatively, the 
Cornwall- Jackiw-Tomboulis (CJT) formalism [31] may be used, as a translationally invariant 
QGP in equilibrium and without spontaneously broken symmetries is considered. In this case, 
both formalisms are equivalent [32] . 

The CJT formalism requires the stationarity of the effective action 

r[D, S] = I - ~ {Tr [inD- 1 ] + Tr [D^D - l] } 

+ {TrflnS- 1 ] +Tr[S c 7 1 S-l]} + r 2 [D,S], (1) 

where / is the classical action and D and S are the full gluon and quark propagators while the 
subscript denotes the respective free equivalents. The functional i~2 represents the sum over 
all two-particle irreducible skeleton graphs of the theory. The traces Tr contain the integration 
over the four-dimensional phase space as well as a trace tr over discrete indices. The integration 
is performed using the imaginary time formalism 33 34 35] . For the grand canonical potential 
n = -TF [36137] this yields 

V = ^/(^M^Mlni}- 1 - II D) +2tr|^n F Mlm(ln S- 1 - SS) |r 2 , (2) 

where tib = (exp (J3u>) — 1) _1 with f3 = 1/Tis the Bose-Einstein, and rip = (exp(/3(w — /i)) + l) _1 
the Fcrmi-Dirac distribution function. 



2.2 Application to QCD 

In our present approach, the infinite sum rjj is truncated at 2-loop order leaving the contribu- 
tions exhibited e.g. in equation (25) in [23]. Within the CJT formalism, the self-energies then 
follow from a functional derivative of J^, giving the well-known 1-loop self-energies (e.g. equa- 
tions (26) and (27) ibid.). In order to achieve a gauge invariant formulation of the model, we 
apply an additional approximation of hard thermal loops (HTL). Although originally being 
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derived for soft external momenta u>,k ~ gT <C T, HTL results coincide with the complete 
one-loop results on the lightcone [23J38] and thus provide the correct limiting behaviour. 

The resulting HTL self-energies can be found in textbooks. Here we follow the conventions 
of Blaizot et al. [15], where essential features of our model have been worked out, and use 



n ltv = n T (u,k) (4r(k)) 
7o ^ = ^+(w,fc) A+(k) 

with the scalar self-energies 
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where M(T, /i, g 2 ) is the thermal fermion mass or plasma frequency and mu(T, /i, g 2 ) denotes 
the Debye screening mass 
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The number of colors N c is fixed at 3. The numbers of light quarks N q = 2 and one strange 
quark, N s = 1, are chosen as in the lattice calculations [28] . 



2.3 Properties of HTL self-energies and dispersion relations 

The real and the imaginary parts of the HTL self-energies (0 and ([6|) are [33] 
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where e(k) is the sign function. The real parts are symmetric with respect to uj, while the 
imaginary parts are antisymmetric and differ from zero only for \u>\ < k, i.e. below the light 
cone. The real and imaginary parts for k — 0.5T are shown in FigurefT] Analogously, the quark 
self-energies fulfill the parity relations ReE + (— uj) = Re£-(uj) and Im£±(—u>) — — ImZ , _(oj) 
as shown for k = 0.5T in Figure El 

It follows directly from eqs. ([8])- ([TT))) that the HTL self-energies do not account for quasi- 
particle widths since the imaginary parts are zero at the poles of the quasiparticle propagators, 
i.e. above the light cone. The nonzero imaginary parts of the self-energies below the light cone 
are due to Landau damping (LD). LD is a collective effect caused by energy transfer between 
the gauge field and plasma particles with velocities close to the phase velocity ("resonant par- 
ticles"). 

Even though the imaginary parts are formally nonzero only below the light cone, retardation 
leads to an infinitely small contribution even above the light cone, giving a definite sign to the 
self-energies for all energies: £(Im77r(w)) = —s(ui), £(Im/7L(w)) = +e(ui) and s(lmS±(ui)) = 
=Fl. Note that this is not related to Landau damping which is found below the light cone only. 



The HTL propagators follow from Dyson's equations as 
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III, and Sj. 1 = — uj ± (k + S±). On-shcll (quasi)particles satisfy a dispersion relation 
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Fig. 1. The real and imaginary parts of the retarded transverse (left) and longitudinal (right) gluon 
self-energies scaled by the Debye mass squared are shown as functions of the energy ui scaled by the 
momentum k which is fixed at k — 0.5T. 




ui/k oj/k 

Fig. 2. The real and imaginary parts of the retarded quark self-energies for the normal (left) and 
abnormal branch (right) scaled by the plasma frequency squared are shown as functions of the energy 
uj scaled by the momentum k which is fixed at k = 0.5T. 




ui/k ui/k 



Fig. 3. The real parts of the inverse gluon propagators D T L scaled by the Debye screening mass 
squared are shown as functions of the energy u scaled by the momentum k which is fixed at k = 0.5T. 



determined by HeD^ L = and ReS 1 ^ = respectively. It is, therefore, useful to first investigate 
the real part of the inverse retarded HTL propagators. 

Due to symmetry properties of £>t l there is just one positive-energy dispersion relation 
above the light cone: uti and uli, respectively. This means that - up to the sign - transverse 
and longitudinal gluons have the same dispersion relations as their anti(quasi)particle counter- 
parts. The additional tachyonic dispersion relation for longitudinal gluons is related to Landau 
damping. Figure [3] explicitly shows the real parts for fixed momentum k = 0.5T. Both D^ L 

are symmetric with respect to to. The zero of KeD^ 1 determines the dispersion relation u>T,k 
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Fig. 4. The real parts of the inverse quark propagators S ± 1 scaled by the fermionic mass parameter 
squared are shown as functions of the energy u scaled by the momentum k which is fixed at k = 0.5T. 




k/m D k 2 /m 2 D 

Fig. 5. The dispersion relations u>T,h of transverse and WL.fc of longitudinal gluon modes scaled by the 
Debye screening mass are shown as functions of the momentum k scaled by the Debye screening mass 
in linear (left) and quadratic (right) scales. 

for transverse gluons. The zero of ReDjJ 1 above the light cone indicates the dispersion relation 
WL,it of longitudinal gluons, while the tachyonic dispersion relation k (below the light cone) 
is due to Landau damping. 

The inverse quark propagators are not symmetric but, as a consequence of the symmetry 
of the self-energy, satisfy the parity property ReS'^ 1 (— w) = — Rc5Z 1 (o;) (cf. Figured]). Hence, 
quarks are described by the positive energy dispersion relation related to S+, while the dis- 
persion relation of antiquarks is found from the negative energy solution of ReSZ 1 = 0. The 
remaining two dispersion relations represent collective quark excitations: the positive energy 
dispersion relation related to S- describes the plasminos, while the negative energy solution of 
ReS^ 1 = represents antiplasminos. Again, a tachyonic solution appears within the regime of 
Landau damping. 

The evolution of the zeros of the real part of the inverse retarded propagators as a func- 
tion of the momentum k gives the dispersion relations w^-. These dispersion relations can- 
not be expressed as analytic functions iv(k) in closed form. ReD^ 1 ^, k, 77i(w,fc)) = and 
Re5,j~ (u, k, Si(u>, k)) = lead to transcendental equations and have to be solved numerically. 
The results are shown in Figures and [SJ Due to the above parity property quarks and anti- 
quarks obey identical dispersion relations up to the sign, as do plasminos and antiplasminos. 



2.4 2-loop QCD entropy 

Given the explicit form of the HTL self-energies and the respective propagators, we evaluate 
the remaining traces tr in eq. ((2]). Assuming equal masses for u and d quarks and zero chemical 
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0.5 1.0 

k/M 

Fig. 6. The dispersion relations uji.u of quarks (solid black), antiquarks (dashed grey), plasminos (black 
dashes) and antiplasminos (black points) scaled by the fermionic mass parameter are shown as functions 
of the momentum k scaled by the fermionic mass parameter in linear (left) and quadratic (right) scales. 



potential of strange quarks the isospin chemical potential fj,j = (/i M — /i<j)/2 is supposed to 
vanish at zero net charge. Therefore, there is only one independent chemical potential /i = fi q — 
/is/3, where \ib is the baryo-chemical potential. As a consequence the flavor trace gives equal 
contributions of light quarks to the thermodynamic potential (i.e. a factor N q ). The contribution 
of the strange quark flavor is here supposed to be equal to the light quark contribution up to 
a substitution of fi — ► and N q — > N s in the following formulae. 

Taking the trace in Minkowski space, the gluonic part decomposes into three contributions 
for one longitudinal and two (equivalent) transverse polarizations, while the quark contribution 
becomes the sum of the normal and the abnormal quark branch (positive and negative chirality 
over helicity ratio, respectively) when taking the Dirac trace. The remaining traces only give 
overall factors: the color trace (N% — 1) for the gluons and N c for the quarks, and the spin traces 
for quarks an additional 2. Defining the prefactors d g = — l,d q 
and introducing the abbreviation J d „ 
reads 



2N r N q and d* 



2N r N, 



Jd"/c/(27r) 4 the HTL grand canonical potential then 
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d g f n B |2im(ln£> T 1 - D T n T ) + Win (-D^ 1 ) + £>l#l)} (11) 
+2 V d, [ n F |lm(lnS+ 1 - S+S+) + Im(ln (-Si 1 ) - Jr 2 . 

Differentiating the thermodynamic potential with respect to the temperature at constant 
chemical potential gives the entropy. In contrast to the pressure, which is influenced by vacuum 
fluctuations, the entropy is sensitive to thermal excitations and therefore manifestly ultraviolet 
(UV) finite. As such, it is ideally suited to investigate the properties of the QGP |15j . 

Due to the stationarity of the thermodynamic potential with respect to the full propa- 
gators, 5Q/5D — 0, only the derivatives of the statistical distribution functions contribute. 
Using Im(£>T-77T) = KeD^lmllT: + ImDxR-eTTT, the entropy density can be written as s := 
-V- 1 3Q/dT\^ = s s , T + s 9 , L + E q , s (si,+ + Si,-) with 



S 3,T 



= -2d r . 



Jd 4 k 01 1 ' 



Sq,h 



1q/s,± 



dnp(u>) 



(12) 
(13) 
(14) 



each describing the entropy density of one quasiparticle species in the absence of the others. 
An interaction (correlation) entropy density contribution would contain terms of the form 
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Imi>rRe77T and the derivative of r^T with respect to the temperature. However, at 2-loop 
order these terms exactly cancel each other [15]. In fact, this seems to be a generic, topological 
feature [39] which has explicitly been proven for QED [40] and <£ 4 theory [41] too. 

We now focus on the terms Imln(±D^ L ) and Imln(±S7_ ), which can be written as 

/ImZ? -1 \ 

ImlnZ^T 1 = arctan E_ + 77£ (ImD^ 1 )e(-ReD^ 1 ) , (15) 

\ Re_D T / 

/ImZ) -1 \ 

Imln (-D h x ) = arctan ^- - ne(ImD^ 1 )0(+ReD^ x ) . (16) 

\RcD L J 

Similar expressions apply for the two quark propagators: one has to substitute ST 1 for D^ 1 in 
and ST 1 for 1 in (fTB|> . 
From the properties of the imaginary parts of the self-energies discussed above, we find 
epmD^ 1 ^)) = — e(u>) for the gluons and e(haS±(uj)) = —1 for the normal and abnormal 
quark branches. We end up with 

s g , T = +2d g [ ^{TC(w)e(-ReD T 1 ) - arctan + ReL» T ImiT T }, (17) 

s s ,l = - d g ^^ £ (uj)0(+R c Dl 1 ) - arctan |^ + Rc£ L Imi7 L }, (18) 

s g/St ± = ±2d q/ „j ^{ nO^ReS^) - arctan |^|| + ReS±Imi; ± }. (19) 

The partial entropy densities (|17p ~(|19 p and, therefore the whole entropy density expression, are 
independent of possible renormalization factors. As required, the expression is also explicitly 
UV finite, as the derivatives of the distribution functions soften the UV behavior. The terms 
7r<9(. . .) represent the quasiparticle contributions to the entropy, while the terms containing the 
imaginary parts of the self-energies arc related to damping effects and quasiparticle widths. In 
the case of HTL self-energies, Landau damping is contained within the latter terms. 

The quark entropy density s q / s = s q / s + + s q / s _ can be simplified by utilizing the parity 
properties for quark propagators and self-energies. Introducing the distribution function of 
antiparticles rip = (e^ u)+l ^ + 1) _1 with dnp(— lu)/8lu = — dn^ {u>)/du> and substituting uo — > — u> 
within s 9) _, we find 

s q/s = 2d q/s (^ + ^) {^(-Rc^ 1 ) - arctan(|||) + ReS + Imi; + } . (20) 

Regarding the quasiparticle pole term nO^ReS^ 1 ), the energy integration from — oo to gives 
the (anti)plasmino contribution, while the integration from to +oo delivers the contributions 
of the (anti)particles to the entropy density. Isolating both parts of the spectrum by applying 
the parity properties once more gives the explicit expressions 

oo 

V»,tl= 2d q/s j ^ (.) |^0(-Re5+ 1 ) - arctan (^§r) + RcS+Im^+J , (21) 



OO 

s q / s ,Pi=-2d q/s J ^ j ^ (.) |7T0( ReST 1 ) - arctan f|^pr) + ReS_Imi;_ 



(22) 



where the sum of the derivatives of the distribution functions is abbreviated by the parentheses 
(.). While this separation seems straightforward, it has to be handled with care as the Landau 
damping term within the quark self-energies £± (see the imaginary parts in Figure [2]) can, in 
general, not be separated into quark and plasmino contributions in this simple way. 
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2.5 The full HTL QPM 

Since the entropy density of the quark-gluon plasma for 2-loop QCD is the sum of the single 
quasiparticle entropy density contributions, it can be considered as mixture of non-interacting 
ideal quasiparticle gases. It is natural to assume that the pressure, which follows from the 
entropy density by integration, consists of single partial pressures, too. Therefore, we use the 
ansatz p = p g ^ + Pg.L + J2i=q s Pi~ B(n^, 27l, E±) for the pressure, where B is chosen appro- 
priately to ensure thermodynamic consistency. The ansatz has to satisfy Sj = dpi/dT\^ which 
leads to 

Pg T = +2d g [ n B (7re(w)6>(-ReL>T 1 ) ~ arctan + ReD T Im7T T }, (23) 

Jd^k RcD- T J 

Pg.L = - d g f n B (7re(w)6>(+Re J DL 1 ) - arctan — — + ReD L Imi7 L }, (24) 
J d 4 k L RcD h J 

p q/s = 2d q/s / (n F +n£)\ 7r6»(-ReS'+ 1 ) - arctan -± + ReS+ImZ' + \, (25) 
Jd 4 k L Reb + J 

where the integrability condition dB/dlli = dp/dlli has to be fulfilled for every quasiparticle 
species i. Thus B ensures the stationarity of the thermodynamic potential under functional 
variation with respect to the self-energies [32] • Note that the plasma frequency within the 
s-quark pressure differs from the plasma frequency within p q as /i s = 0. 

The pressure fully defines the model. The particle density follows by differentiation of the 
pressure with respect to the chemical potential at constant temperature. The Bose-Einstein 
distribution function «b does not depend on [i g and strange quarks are included into the model 
with manifest zero net particle density, therefore n g ^ = = n s = . Due to the integrability 
condition, the terms containing the derivatives of the self-energies with respect to fi vanish, so 
that ^ 

n=n q =2dJ f^ + ^F ^Lot-ReST 1 ) -arctan +Rc5+Imr + 1 , (26) 

id*fe\ oil dfi J I ReS~ + L J 

thus n q (n — > 0) — »■ 0. 



2.6 Effective coupling 

Obviously, 2-loop QCD is only a crude approximation of the full theory. In order to ac- 
commodate non-perturbative effects in the quasiparticle model, we introduce some flexibil- 
ity by parameterizing the QCD coupling constant g 2 in a phenomenologically motivated way. 
The truncated 2-loop running QCD coupling g 2 is given by g 2 {x) — 167r 2 (/3 ln(x)) -1 (I — 
2f3 x ln[m(»] (^ln(x))- 1 ) 13). where f3 = llN c /3 - 2JV//3, 0i = 51 - l9N f /3 and N f = 
N q + N s . It depends on the ratio x = jl 2 / A 2 of the renormalization scale p, and the QCD scale 
parameter A. A term involving ln _2 (/2 2 /yl 2 ) which is only a small correction for p 2 ks A 2 was 
neglected. 

The renormalization scale is usually taken to be the first Matsubara frequency 27rT, while 
the latter one is just a parameter to be adjusted using experimental data. Introducing the 
pseudocritical temperature of QCD matter at vanishing net baryon density T c and substituting 
A — » 2ttT c /\, the ratio fl/A becomes XT/T C . In order to avoid the Landau pole of g 2 (T/T c ) at 
T c a temperature shift with parameter T s is introduced, replacing x by £ = A(T — T s )/T c . The 
result 

G 2 (T >T u-0)- J^f_ ( i ^ihMl) (27) 

is our effective coupling. Within the plasma phase for temperatures T > T c /A + T s it is well- 
behaved; however, at some point within the hadronic phase, i.e. below T c , an infrared (IR) 
divergence does occur. In order to prevent this divergence a phenomenological infrared cutoff 
for G 2 can be applied. See [44] for details. 
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2.7 Adjustment to lattice calculations 

The two QPM parameters A and T s have to be adjusted to results of numerical first-principle 
QCD calculations dubbed lattice data. Most of the past work on the QPM has been tested 
against lattice data from [45] for rather large and temperature dependent lattice restmasses 
of m q — OAT and m s = LOT compared to the physical quark masses m Uj( i ~ lOMeV and 
m s ~ 90 — 150 MeV [43] . Recently, new lattice data has become available [28], which relies on 
lattice restmasses much closer to the physical quark masses and which is used in this work. 

Also, lattice calculations are performed on a finite lattice, while our quasiparticle model 
is formulated in the thermodynamic limit, i.e. aimed at describing a spatially infinite plasma. 
In order to compare our model with lattice data, the proper continuum extrapolation of the 
latter one is required. A safe continuum extrapolation on the lattice is a fairly demanding 
work. Therefore, various estimates have been applied, e.g. simply scaling the lattice results by a 
factor being strictly valid only for asymptotically high temperatures or for the non-interacting 
limit. To account for a possible deficit of such rough continuum estimates of the lattice data 
we introduce an ad hoc scaling factor d\ a t which turns out to be nearly unity. 



2.8 Nonzero chemical potential 

The parametrization of G 2 in eq. ([27]) is valid for \i = only. However, it is possible to use the 
thermodynamic consistency of the QPM to map the results at zero chemical potential into the 
T-fj, plane [46]. Specifically, this means to impose the Maxwell relation ds/dfi\T = dn/dT\^ on 
the thermodynamic quantities. Ordering the terms with respect to the partial derivatives of the 
effective coupling gives an elliptic quasilinear partial differential equation 

dG 2 dG 2 , 
aT llT +a ^ = b > (28) 

named hereafter flow equation, with the coefficients ax, and b depending on T and fi explicitly 
and, via the self-energies, implicitly. It is solved by the method of characteristics by introducing 
a curve parameter x, assuming that T = T(x), fx = (i(x) and G 2 — G 2 (x). Subsequently, the 
comparison of G 2 X = G 2 T T :X + G 2 ^n. x with the flow equation gives a system of three linear, 
coupled ordinary differential equations: G\ — —6, T rX = —cit and [i. x = — &n which can be 
solved using standard numerical methods. The initial condition for the flow equation is the 
effective coupling at /i = 0, with model parameters fixed by comparison of the entropy density 
with lattice results. 



3 The effective QPM 

Assuming that transversal gluons and quark particle excitations propagate predominantly on 
mass shells the full HTL QPM can be significantly simplified as it implies explicit (asymptotic) 
dispersion relations Wj(fc) of the form uij(k) = k 2 +m| QO as approximations to the full, implicit 
ones [18] . The mi j00 terms depend neither on energy nor momentum and are therefore called 
asymptotic (thermal) masses. In order to adjust the eQP to lattice data, they can be modified to 
accomodate lattice restmasses Wi using a prescription from [47]: rr? ioo — > m 2 +2mimi i00 +2rn 2 oa . 
Additionally neglecting Landau damping, i.e. assuming vanishing imaginary parts of the self- 
energies, and contributions from collective excitations, i.e. plasmons and (anti)plasminos, leads 
to the eQP. As collective excitations are exponentially suppressed^! and the effect of Landau 
damping is small at vanishing chemical potential the eQP seemed to be sufficient. 

That is after calculating the propagators using Dyson's relation from the HTL self-energies, the 
residues of the poles in the spectral density of both plasmon and (anti)plasmino propagators vanish 
exponentially for momenta k ~ T, [i, which give the dominant main contributions to thermodynamic 
integrals. 
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Fig. 7. The scaled entropy densities s/T 3 of the full HTL QPM with quadratic IR regulator (solid 
black lines; T s = 0.728T C and A = 6.10) and the eQP (grey dashed lines; T 3 = 0.752T C and A = 6.26) 
adjusted to lattice data for Nf = 2+1 from [28] with di at = 0.96 are shown as functions of the scaled 
temperature T/T c . The adjustment quality of the full HTL QPM to lattice data is indistinguishable 
from the eQP. The single contributions to s HTL , including their respective LD contributions, are given 
in the right figure (dashed black: transversal gluons+(anti)quarks, dash-dotted: longitudinal gluons, 
dotted grey: (anti)plasminos). 



However, as previous studies of the flow equation 18 44J have shown the characteristic 
curves emerging at T rj T c cross each other in some region of finite values of /i for parameters 
adjusted to lattice QCD results (cf. dashed lines in Figure [5] below) . This unfortunate feature 
prevents an unambiguous extrapolation of thermodynamic quantities into the full T-^-plane. 
Romatschke [27] has shown for the 2 flavor case that these crossings can be avoided by using 
the full HTL model. We are going to extend his line of arguing to the physically interesting 
case of 2 + 1 flavors. 



4 Investigation of the full HTL QPM 

4.1 Adjustment to lattice data at fi = 

While the eQP is able to accommodate arbitrary lattice restmasses by means of a modified 
asymptotic dispersion relation, the full HTL model relies on the HTL dispersion relations 
and thus massless particles. We assume here that the employed lattice restmasses in [28] are 
sufficiently small to be absorbed in suitably adjusted parameters. 

For T > T c and Nf = 2 the full HTL model has been shown to give a description of lattice 
data being equally well as the eQP [48] . For Nf = 2 + 1 flavors we meet a similar situation: 
The full HTL QPM describes the lattice QCD data [28] as good as the eQP model (see Figure 
left). The extension to T < T c , on the other hand, is not straightforward. Instead of a linear 
IR regulator [49] , it is necessary to use a quadratic parametrization of the effective coupling in 
order to achieve agreement with lattice data also below the pseudocritical temperature, which 
we consider here as mere parametrisation of the lattice data. 

When evaluating the individual contributions to the entropy density of the full HTL QPM we 
find the entropy density contributions of longitudinal gluon s 9i l (eq. (I18p ) and (anti)plasminos 
s q,p\ ( cc f ([22]) ) to be negative. This is due to the fact that both represent collective phenomena of 
the QGP resulting in correlations not present in a noninteracting medium. As a consequence, the 
transverse gluon and (anti)quark entropy density contributions have to increase in comparison 
to the eQP in order for the sum of the partial entropy densities to describe the same lattice 
data as the eQP (see Figure right). This not only allows for a pure quasiparticle entropy 
contribution much closer to the Stefan-Boltzmann limit than in the eQP but also proves to 
have a positive impact on the extension to nonzero chemical potential. 
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Fig. 8. The solid curves in both graphs are several characteristics of the full HTL flow equation for 
2 + 1 quark flavors using parameters from the adjustment of the full HTL QPM to lattice data from 
[28j shown in Figure [7] The characteristic curve emerging from T c is depicted as bold solid line. All 
crossings have disappeared (left panel). For a comparison, the right panel shows the characteristics 
of the eQP flow equation using the parameters of the eQP adjustment to the same lattice QCD data 
(dashed curves). 



4.2 Solution of the flow equation 

Solving the flow equation (J2HJ) with coefficients listed in Appendix [A] the characteristics are 
found to be well-behaving, as can be seen in Figure [5) Also a stronger curvature of the full HTL 
characteristics compared to the eQP characteristics (shown as dashed lines in the right panel) 
is observable. 

To explain the disappearance of the ambiguities caused by crossing characteristics in the 
eQP we mention that the crossings appear due to the effective coupling G 2 being too large 
near the pseudocritical temperature [50]. Since the entropy density increases with decreasing 
mass parameters m 2 D and M 2 (which are proportional to G 2 T 2 at /.i = 0) the crossings would 
therefore disappear for a larger eQP entropy density. One way to allow for a larger eQP entropy 
density is to take into account collective modes. As medium effects indicate correlations between 
the constituents of the eQP plasma, including them causes a decrease of overall entropy density. 
Consequently, the eQP parameters have to change in order to still describe the same lattice data, 
causing the entropy density to increase. With the resulting decrease of the effective coupling 
G 2 the crossings partially disappear. 

However, the different parametrization at fx = alone cannot account for the complete 
absence of ambiguities for the full HTL model. Instead, the influence of collective modes and 
Landau damping on the flow equation has to be examined. We therefore calculate the character- 
istics of the full HTL model respectively disregarding terms stemming from these contributions. 
While neglecting plasmon/(anti)plasminos terms from the coefficients ax, a M and b (see Ap- 
pendix [X| but keeping the Landau damping contributions leads to deformed characteristics 
meeting T — at smaller /i and no crossings appear. Hence, it is neither the plasmon nor 
the plasmino term which accounts for the vanishing crossings. However, neglecting the Landau 
damping terms immediately leads to crossing characteristics. Therefore, both collective exci- 
tations (in order to obtain a reasonably small coupling G 2 ) and Landau damping (in order to 
ultimately remove the crossings) are necessary to obtain a flow equation with unique solutions. 
Using this flow equation, it is possible to extrapolate the equation of state from lattice QCD 
at fi = towards T = 0. 



5 Conclusion 

The mapping of a previous quasiparticle model (eQP) into the T-fi plane was plagued by 
crossing characteristics. It is shown here for the 2+1 flavor case that, if using the full HTL 
model, these crossings disappear. Collective modes (longitudinal gluon and (anti)quark hole 
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excitations, i.e. plasmons and (anti)plasminos respectively) as well as Landau damping of the 
collisionless quasiparticle plasma, both neglected hitherto in the eQP, need to be taken into 
account to avoid the ambiguities. 

With the problem of crossing characteristics solved, one can proceed to derive an equation of 
state, following from the full HTL QPM, especially for the cold and dense quark-gluon plasma 
of interest in future heavy ion collision experiments or for the simulation of possible quark stars. 

Acknowledgment: R.S. would like to thank the organizers of the Zimdnyi 75 Memorial Workshop for 
the invitation to present his results at this very inspiring workshop. 



A Coefficients of the flow equation 

For the reader's convenience and to extend the results in [48 to Nf — 2 + 1 flavors the full 
HTL flow equation is presented. To calculate the Maxwell relation, the derivatives ds g /d/j = 
A g dm 2 D /dii, (ds q /du) impL = A q dM 2 jdix, (ds s /du) impL = A s dM 2 /dii and (dn q /dT) impL = 
A n dM 2 J dT are necessary. The explicit derivatives cancel within the Maxwell relation due to 
Schwarz's Theorem. We find 



2ir 3 m 2 D 



dkk 2 



do> 



dn B 4(w 2 - k 2 )lm 3 n T 
~df~ (Re 2 D- 1 +Im 2 7T T ) ; 



dni 
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with abbreviations At := d/dT(np + rip) and := d/du(rip + rip). The derivative of the 
strange quark entropy density with respect to the temperature equals the light quark expression 
at vanishing chemical potential with A s — A q (fi — 0) and dM 2 /da = (<9M 2 /<9/i)| M= o. Imposing 
the Maxwell relation ds/dfi\r = dn/dT\n and employing the prefactors in eq. ((7|) the coefficients 
of the flow equation ([28]) are given by 



N 2 



1 



lWr. 



It- 



Am 



1 A n 9 N c N q 9 

- 2A^ C + Nq + N s ] T + -ArV 

6 V TT Z 



+- 



A^ 2 



1 



16JV C 



A n 



m-i . 



~-T 2 A s 



b = 



N 2 



1 



8JV r . 



TG 2 A r , 



Nr.N, 



3tt 



N 2 



8N r TT 



\aG 2 A q . 



(32) 
(33) 

(34) 
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These expressions correct a few typos in [48] (see eqs. (B.1)-(B.5) therein). 
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